RPA Analysis of a Two-orbital Model for the BiS2-based Superconductors 
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The random-phase approximation (RPA) is here applied to a two-orbital model for the BiS2-based supercon- 
ductors that was recently proposed by Usui et al., arXiv: 1207.3888 Varying the density of doped electrons per 
Bi site, n, in the range 0.46 < n < 1.0, the spin fluctuations promote competing A\g and B2g superconducting 
states with similar pairing strengths, in analogy with the A\g-Big near degeneracy found also within RPA in 
models for pnictides. At these band fillings, two hole-pockets centered at (0, 0) and (tt, tt) display nearly paral- 
lel Fermi Surface segments close to wavevector (7r/2, 7r/2), whose distance increases with n. After introducing 
electronic interactions treated in the RPA, the inter-pocket nesting of these segments leads to pair scattering with 
a rather "local" character in k-space. The similarity between the A\g and B29 channels observed here should 
manifest in experiments on BiS2 -based superconductors if the pairing is caused by spin fluctuations. 

PACS numbers: 74.20.Mn,74.20.Rp,74.70.-b 



Introduction. — The recently discovered family of layered 
bismuth oxy-sulfide superconductors^^ has immediately at- 
tracted considerable attention from the Condensed Matter 
community due to its close similarities with the famous iron- 
pnictide superconductors;^— As in the case of other layered 
unconventional superconductors, such as the cuprates and the 
aforementioned iron pnictides/chalcogenides, this new fam- 
ily displays a layered structure involving BiS2 planes where 
the observed superconductivity is believed to reside. The 
first report of superconductivity originated in 6140483, with 
Tc = 4.5 Kji Superconductivity has also been reported in 
i?eOi_2:Fj:BiS2, where Re = La, Nd, Ce, and Pr, with cor- 
responding Tc = 10. 6^^ 5.6^^ 3.0r^ and 5.5 KJ^ These com- 
pounds are metallic in the normal state and Density Functional 
Theory calculations indicate that the relevant bands crossing 
the Fermi surface (FS) originate mainly from the Bi 6p or- 
bitals, as shown, e.g., for LaOi_a;Fa;BiS2'^ However, con- 
trary to the majority of the Cu- and Fe-based unconventional 
superconductors, no magnetically ordered phase has been de- 
tected thus far in the BiS2 compounds. This apparent ab- 
sence of magnetism in the BiS2 compounds may still locate 
them in the same category as LiFcAs, FeSe, and possibly 
Sr2V03FeAsrl that are also non magnetic but their pairing 
properties are widely believed to still originate in short-range 
magnetic fluctuations. For these reasons, and despite the ab- 
sence of observed long-range magnetism in BiS2, it is impor- 
tant to study the potential role of spin fluctuations in these 
novel materials and the pairing channels that those fluctua- 
tions tend to favor, to help in the analysis of experimental data. 

In this manuscript, the two-orbital (2-orbital) model re- 
cently introduced by Usui et al. is adopted^ The fact that 
the relevant orbitals in BiS2 compounds are p-type, where 
Coulomb interactions should be smaller than in d orbitals, 
turns RPA into a suitable technique, whose results deserve a 
careful analysis if electron correlations are found to be impor- 
tant for superconductivity in these materials. Similar calcula- 
tions for a related four-orbital model"* are underway. Note that 
in Ref. H a brief discussion of RPA calculations has already 
been presented. The results discussed by Usui et al. consisted 



of a single set of couplings (equivalent to our J/U = 0.2 cal- 
culations below) at n = 0.5. Their early weak-coupling RPA 
analysis is here expanded via a systematic study of the influ- 
ence of the band filling n and the identification of the domi- 
nant channels for superconductivity under the assumption of 
a spin fluctuations mechanism. The main novel contribution 
of our present effort is the identification of closely competing 
B2g and Aig gap functions as the dominant pairing channels, 
particularly for band fillings around n = 0.5. At quarter fill- 
ing {n = 1.0), another pair of almost degenerate gap func- 
tions (with symmetries and Big) is found to closely com- 
pete with the previously mentioned dominant pair, especially 
at J/U = 0.3. 

Hamiltonian. The 2-orbital model described by Usui et ali^ 
contains hopping parameters up to fourth neighbors, and in 
k-space is given by 
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4^2^^ (cos 2kx — cos 2ky 



+ 4:t2x+y (cos 2kx cos ky — cos kx cos 2ky) . 

The operator dl^ ^ ^ {du.L>.cr) in Eq. ^ creates (annihilates) an 
electron in band f = X,Y, with spin ct = ±, and wavevec- 
tor k. The values for the hopping parameters are those from 
Ref. m and are reproduced in Table I for completeness (in eV 
units, as used throughout this paper). Figure [Tt a) shows the 
FS hole-pockets for four different band fillings n = 0.46, 
0.5, 0.65, and 1.0, with corresponding chemical potentials 
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TABLE I: Tight-binding parameters (eV) for 2-orbital model. 
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FIG. 1 : (Color online) (a) Hole-pockets for four different electronic 
fillings: n — 0.46 (solid red), n — 0.50 (dashed green), n = 0.65 
(dotted blue), and n = 1.00 (dot-dashed magenta). Note that close 
to the (7r/2,7r/2) wavevector, where the n — 0.46 pockets almost 
touch, the increase of n decreases the radius of the hole-pockets and, 
more importantly, the adjacent FS segments (inside the dashed box) 
become more and more parallel, (b) Lindhard function xo for the 
same fillings as in panel (a). Note that the position in k-space of the 
leftmost peak is clearly associated to FS nesting through a (fc„, 0) 
vector, as indicated in the inset, which zooms-in the dashed box in 
panel (a). Indeed, the position of the leftmost peaks in xo agree 
(within a few percent) with the vectors indicated in the inset (see text 
for details, especially Fig.O. Obviously, there are additional nesting 
vectors that become evident in a 2-d plot of xo [Fig.|5jb)]. 



fi = 1.10375, 1.12514, 1.21828, and 1.52621 (in principle, 
71 = a; in LaOi_a;Fa;BiS2).'^ Panel (b) shows the correspond- 
ing non-interacting magnetic susceptibilities xo- The leftmost 
peaks in xo^ located at (fc„, 0), with < fc„ < tt/2 as the fill- 
ing varies from n ~ 0.46 to 1.0, can be associated to FS nest- 
ing once it is noticed that their position matches the horizon- 
tal separation between the two adjacent FS segments from the 
pockets centered at (0, 0) (F) and (tt, tt) (A/), as highlighted 
by the dashed box in panel (a) and sketched in the inset to 
panel (b). Note that the horizontal separation is well defined if 
the two FS segments are parallel, which is the limiting case as 
n increases, as shown in the inset, to n = 1.0 (for details, see 
Fig.|5]and the associated discussion). It is also important to re- 
mark that once interactions are introduced, the leftmost peak 
in Xo is the one that diverges in the RPA calculation of the 
spin susceptibility xrpa for almost all the fillings and various 
values of interaction parameters. This divergence indicates a 
tendency to magnetic order, or at least strong spin fluctuations 
(paramagnons), with characteristic wavelength determined by 
(fc„,0). Our analysis is not extended into the n < 0.45 re- 
gion since there the topology of the FS changes (see Ref. [U 
for details of the FS at lower fillingsa^). 



FIG. 2: (Color online) RPA spin susceptibility (solid red curves in 
the main panels) and dominant gap function (red and blue dots in 
the insets) for (a) n = 0.46 and (b) n = 0.50. In the inset to each 
panel, the dominant gap function with symmetry B29 is shown. The 
subdominant gap function (not shown) has symmetry Aig and its 
eigenvalue is almost degenerate with the dominant one (see text). 



The Coulomb interaction in the Hamiltonian is given by 

Hint = U ^ni,a,^ni^a.i + [W - J/2) ^ rii^anifi 

i,a i,a</3 
i,a<l3 

+ •^11 (4«,t4a,4'^i,/3,i'^i,/3,t + ^-C-), 

i,a<f3 

where the notation is standard and the many terms have been 
described elsewhereiS Here, the usual relation U' = U~2J 
is assumed, and J/U is a parameter. Calculations were done 
for 0.1 < J/U < 0.4, in steps of 0.1, for the four fillings n = 
0.46, 0.50, 0.65, and 1.0. The multi-orbital RPA calculations 
performed here follow closely those described in Ref.H^ and 
previous works by the authorSf^i^ All results were obtained 
at temperature T = 10^^ and an imaginary part 77 = 10~^ 
was used to regularize the Green's functions. 

Our RPA results for spin-singlet pairing link the dominant 
superconducting gap functions to spin fluctuations, which 
originate in FS nesting and are enhanced by electronic inter- 
actions. The particular relative topology of the two adjacent 
hole-pockets (see Fig. [Hi promotes pairing whose strength is 
independent of the global symmetry of the pairing functions 
[see Fig. Efb)]. Indeed, the i?2g and Aig symmetries have 
essentially the same pairing strength, which is determined 
by pair scattering between these two adjacent FS segments 
(see Fig. Is) close to {tt/2, tt/2) in the Brillouin Zone (BZ). 
In addition, our results show that both dominant gap func- 
tions change sign between these two segments (Figs. |2]to|4l), 
and the pairing is through the intraorbital scattering channel 
[Fig.[3jb)]. The near degeneracy Aig-B2g is the analog of the 
near degeneracy Aig-Big found also in RPA calculations for 
the pnictidesf - since the pocket structures in both cases can 
be related by a 45° rotation. Results for spin-triplet pairing 
are presented in the supplemental material at the end of the 
manuscript. 

Results and Discussion. As mentioned above, the most im- 
portant feature of the FS for fillings between 0.46 and 1.0 is 
that the hole-pockets centered at the T and AI points present 
almost parallel segments close to the {tt/2, tt/2) wavevector. 
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FIG. 3: (Color online) (a) RPA spin susceptibility and dominant gap 
function for n = 0.65. Orbital composition for the (0, 0) and (tt, tt) 
FS pockets = 0.65), (c) and (d), respectively. The winding angle 
Q is counter-clockwise, starting from the k^^ direction. Assuming the 
nesting described in the inset to Fig.[TJb) as producing the spin fluc- 
tuations that provide pairing, the pair coupling is then intraorbital. 
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FIG. 4: (Color online) (a) Dominant gap function with symmetry 
A\g at n = 1.0. (b) Main panel: normalized pairing strengths A for 
the dominant {A\g, solid red curve) and subdominant (B2g, dashed 
green curve) gap functions. Although the two curves are very close, 
the eigenvalues are not degenerate. In the inset, the structure of the 
subdominant gap function {B-zg) is shown. When compared to that 
of the dominant one \_A\g in panel (a)], it is clear that the structure 
around (7r/2, 7r/2) is very similar for both of them, explaining why 
the pairing strengths (eigenvalues) are the same. The region inside 
the dashed box, in panel (a), is analyzed in detail in Fig.[5ja). 



becoming more and more parallel as the pockets shrink, with 
increasing filling [see Fig. [Ha) and inset in Fig. [Tfb)]. In 
Figs. |2] and [3] it will be shown that this has important con- 
sequences for the spin fluctuations and the superconducting 
pairing associated to this 2-orbital model. Indeed, as dis- 
played in the main panel of Fig. [2a) (solid (red) curve), there 
is a divergence in the RPA spin susceptibility for very small k 
values: ko.46 ^ (7i'/25, 0) for n = 0.46, and ko.50 ^ (i'/8, 0) 
for n = 0.50 [panel (b)]. A divergence in the spin susceptibil- 
ity xrpa may point to magnetic order, or at least to strong spin 
fluctuations with wave vector k„. Figure[3ja) shows the same 
calculations, but now for n = 0.65. Note that although xo 
displays a broad-peak structure around (tt, tt) [see Fig. [lib)], 
Xrpa does not present a divergence in this region. In the in- 
sets to Figs. [2ta) and (b), and Fig. [3ja), it is shown that the 
dominant gap function at the FS has symmetry for the 
three cases, showing that despite the changes in the size of 
the hole-pockets the results are qualitatively the same. Fig- 
ures[3jb) and (c) contain the orbital contribution {X, red solid 
curve; Y , green dashed curve) of the BZ states at the FS for the 
r and M pockets, respectively. It is interesting to note that the 
modifications in the position of the peak in xrpa correlates 
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FIG. 5: (Color online) (a) Region around point (7r/2, 7r/2) of the BZ 
[dashed box in Fig.Oa)], showing the dominant (Aig) gap-function 
for n = 1.0 and J jlj = 0.3. (b) Two dimensional contour plot of 
Xo also for n = 1.0. The horizontal (blue) vector in panel (a) con- 
nects the maximum amplitude of the gap-function in both pockets. 
Note also the horizontal (blue) vector in panel (b) along the k^^ direc- 
tion, indicating the position of the maximum value of xo- These two 
vectors agree up to a difference smaller than the width of this maxi- 
mum peak in xo- Therefore, it can be shown (see text) that the line 
describing the position of the points in the M pocket in relation to 
the points in the V pocket, as indicated by the two additional vectors 
(black and red) in panel (a), satisfies ky ~ — fc^ -f where (fc„ , 0) 
and (0, kn) are the positions of the maxima in xo (with n — 1.0). 
This equation also describes the line of local maxima of xo, as seen 
in panel (b), originating from FS nesting. 



well with the "separation" between the T and AI hole-pockets 
in the region around (7r/2, 7r/2). For the purposes of describ- 
ing our results, this separation will be defined as the horizontal 
distance between two parallel lines tangent to the hole-pockets 
at the points where each intercepts the F — M (S) line. As de- 
scribed in more detail in Fig[5la) [and already mentioned in 
connection with Fig[TJa)], as the filling increases these seg- 
ments of FS approach more and more the parallel lines just 
defined, justifying the definition just given. 

The RPA results for the gap functions also point to an inter- 
esting effect, namely, the small value of k„ for fillings 0.46 < 
n < 1.0 results in the pairing strength depending on very "lo- 
cal" properties of the gap function at the adjacent segments of 
the hole-pockets. This implies that the pairing strength of gap 
functions with different symmetries is very similar, as long as 
they have the same "local" properties. To demonstrate that, 
in Fig. [Ifa) the dominant gap function (with Aig symmetry) 
is shown for n ~ 1.0 and J/U ~ 0.3. It is clear that this is 
very similar in structure to the subdominant one shown in the 
previous figures. In the inset to Fig.[4tb) the subdominant gap 
function with symmetry B2g is displayed for the same param- 
eters. Comparing it with the dominant gap function in panel 
(a) note that, despite having different symmetries, the two gap 
functions are identical in the two adjacent hole-pocket seg- 
ments that cross the E line. For this reason, their pairing 
strengths as measured by A (the eigenvalues of the Eliashberg 
Equation), and shown in the main panel of Fig. [Hb), are the 
same to the third decimal place. Note that the two eigenvalues 
for symmetries Aig and B2g are not degenerate. This seems 
a strong indication that the "local" aspect of the pair scatter- 
ing, as mentioned above, seems to be determinant to estab- 
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FIG. 6: (Color online) Two-dimensional plot of the RPA spin suscep- 
tibility for n — 1.0. The parameter values are J/U = 0.3 and U — 
1.65. The similarity to the results in Fig.|5tb) is clear, showing also 
that there are relevant nesting features along the ky — —k^ + ki.o 
line. Note that a smaller value of U than in Fig. |5la) was used to 
avoid having a peak at (fci.o, 0) that would wash out the features in 
the rest of the BZ. 

lish the pairing properties of this model, at least in our RPA 
weak-coupling approach. It should be noted that the eigen- 
value results shown in Fig. |4]are basically identical to those 
for lower fillings, shown in previous figures, with the only 
difference being the order of the dominant and subdominant 
symmetries. Since their eigenvalues are almost identical for 
all fillings studied, this does not have a special significance. 
Note that xrpa for n = 1.0 and J/U = 0.3 (not shown) 
follows the same trends as described in Figs. |2] and [3] From 
the orbital composition in Fig. Ob) and the gap structure in 
Fig- Hit appears that the symmetry of the and Ai^ pairing 
operators is determined by the orbitals, while the spatial form 
in both cases is characterized by symmetric nearest-neighbor 
pairing with rotational invariance. Thus, the pairing operators 

have the form At = /(k)(4,x,t^^-k,X4 ± dl rtdl^^.Y^) 
where the + (— ) sign corresponds to Ai^ (62^) symmetry 
with /(k) — cosfcj; + cosfcj,, plus higher harmonics with 
Aig symmetry. 

Figure|5|a) shows in more detail the almost parallel FS seg- 
ments of the two hole-pockets for n = 1.0. In this figure, the 
horizontal (blue) vector that was defined above as the separa- 
tion between the two FS segments is displayed. A vector with 
the same length is reproduced in panel (b), where a 2d plot 
of xo in the first quadrant of the BZ is also shown. It clearly 
indicates that the position k„ of the main peak in xo is exactly 
given by the horizontal separation. Not only that, the (red) 
vector along the E line in panel (a) is also reproduced in panel 
(b) and it coincides also exactly with a local maximum of xo- 
In fact (see in both panels the black vectors located at angle 
a), the locus of the ridge of local maxima in xo in panel (b) 
exactly coincides with the BZ points defined by the vectors 
connecting the two FS segments for < a < 7r/2. Figure |6] 
shows the RPA spin susceptibility for n = 1.0. The similarity 
between these results and those in Fig. |5ja) is clear, indicat- 
ing that the FS nesting for the interacting system is the one 



described by the vectors in Fig. |5] Finally, an important issue 
should be highlighted; the four points in the hole-pockets in 
Fig.|5ja) where the gap function has a very pronounced peak, 
are exactly the two pairs of points (one in each pocket) con- 
nected by (fci.0,0) and (0, fci.o). This fact clearly links the 
pairing properties with the spin fluctuations. Note also that 
for n = 1.0 and J/U = 0.3, the second pair of eigenvalues 
(A3 = 0.9038 and A4 ~ 0.9036) corresponds to symmetries 
and Big, respectively (not shown). The same occurs for 
J/U = 0.1 and J/U = 0.2, also for n = 1.0 (but the eigen- 
values are smaller). Yet, the same explanation as described in 
Fig.|5]applies. See the supplemental material for a connection 
between the emergence of a Big symmetry at n = 0.50 with 
the one-dimensionality of the bands. 

Conclusions. Summarizing, a weak-coupling RPA analy- 
sis of a minimal 2-orbital model was used to investigate the 
pairing properties of the BiS2-based superconductors. Fill- 
ings between 0.46 and 1.0 were analyzed. The Hund's cou- 
pling was varied in the range 0.1 < J/U < 0.4. Qualita- 
tively, the results are similar for all values of J/U and differ- 
ent fillings. In the RPA results described here, a clear rela- 
tionship is found between quasi FS nesting, spin fluctuations, 
and superconductivity: the topology of the two hole-pockets 
is such that they present almost parallel segments close to the 
{tt/2, tt/2) wavevector in the BZ. It is found that the horizon- 
tal distance (fc„, 0) between the tangents to these segments at 
the points where they cross the E line is also where the non- 
interacting susceptibility xo has a pronounced peak at (/s„, 0), 
for 0.46 < n < 1.0. Once interactions are introduced, this 
peak will diverge at a certain critical coupling U for each fill- 
ing, and all the values of J/U studied (with exception of one: 
n = 0.5, J/U ~ 0.1). In addition, a line of local maxima, 
connecting the BZ points (fc„,0) and (0, fc„), is clearly ob- 
served in a 2-d plot of xo- As expected, this line can also be 
associated to FS nesting. This nesting structure gives origin to 
pairing functions with similar eigenvalues, i.e., similar pairing 
strengths, and symmetries i?2g and Aig. This close competi- 
tion originates in the FS quasi nesting properties, which deter- 
mine the spin-fluctuation-mediated inter-pocket pair scatter- 
ing. This pair scattering is overwhelmingly between two adja- 
cent FS segments, therefore the properties of the pairing func- 
tions, including the pairing strength, are quite "local", hav- 
ing almost no dependence on their global symmetry. One can 
then predict that pairing symmetry measurements may con- 
tain a mixture of both symmetries if the pairing mechanism is 
driven by spin fluctuations. 

GBM acknowledges fruitful conversations with K. Kuroki, 
Q. Luo, and H. Usui. ED and AM were supported by the 
National Science Foundation Grant No. DMR-1 104386. Af- 
ter finishing this manuscript, a related effort addressing the 
pairing symmetry of these materials using a spin model was 
published^ There, it is found a dominant Aig state analogous 
to ours, but no competing B2g state. 
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FIG. S 1 : (Color online) (a) Dominant and (b) subdominant gap func- 
tions in the spin-triplet channel, for parameters n = 0.5, J/U = 0.3, 
and U = 1.416. 



I. SUPPLEMENTAL MATERIAL 

Spin-triplet pairing. We also tested the two-orbital model 
for the case of spin-triplet pairing. Using the same RPA all the 
four fillings studied in this work were investigated, but calcu- 
lations were carried out only for J/U = 0.2 and 0.3. All 
the critical values obtained for the Hubbard repulsion U were 
slightly above those obtained for the singlet pairing channel. 
However, they were close enough to warrant a brief discussion 
in this supplemental material. Note that in Usui et alJ- the pos- 
sibility of spin-triplet pairing was mentioned, in connection 
with the similarity of the BiS2 bands with those of Sr2Ru04, 
in regards to their common one-dimensionality character Fig- 
ure [ST] shows the gap functions [dominant in panel (a) and 
subdominant in (b)] for parameters n = 0.5, J/U = 0.3, 
and U = 1.416. This critical value of U should be compared 
with that obtained for singlet pairing for the same parameters 
(i.e., U = 1.363, see Fig.|2|b) in the main text). The symme- 
tries for both the dominant (i?2c/) and subdominant (Aig) gap 
functions in the spin-triplet channel are the same as for the 
spin-singlet channel. The main difference here is that they do 
not have as competing pairing strengths as in the spin-singlet 
channel. Indeed, the eigenvalues for Fig. [ST] are Ai = 1.00 
and A2 = 0.88, while for the same parameters in the spin- 
singlet channel their values are Ai 0.989 and A2 = 0.985. 

Quasi one-dimensionality. As mentioned in the main text, 
the two-orbital model has a quasi one-dimensional (Id) char- 
acter, with the hopping between next-nearest-neighbor being 



dominant (tj-j-j, = 0.88, see the Table in the main text con- 
taining the hoppings). It is then interesting to verify how 
the results are modified in case all the other hoppings are re- 
moved from the two-orbital model Hamiltonian, except for 
— 0.05. The energies of the orbitals were kept the same 
as in the original model. RPA calculations for the spin-singlet 
pairing channel were done for n — 0.5 (with correspond- 
ing chemical potential /i = 1.18037278), J/U = 0.2, and 
J/U = 0.3. In addition, the spin-triplet pairing channel was 
investigated for J/U = 0.3, but, again, the critical value ob- 
tained for the Hubbard U was higher than for the singlet chan- 
nel, therefore, these results are not shown. Singlet pairing 
results for both values of the Hund's coupling were similar, 
therefore, just the results for J/U = 0.3 will be presented. 
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FIG. S2: (Color online) Band structure for a quasi Id Hamiltonian 
obtained from the two-orbital model discussed in the main text by 
keeping only two hopping terms: t'^±^ and t^^ (see text for details). 
The Fermi energy is at Ef ~ 0.0. 




FIG. S4: (Color online) (a) Dominant and (b) subdominant gap func- 
tions in the singlet channel for the quasi Id model. The parameters 
used are n — 0.5, J/U = 0.3, and U = 0.8086. The arrows in- 
dicate pair scattering processes associated to the peaks in Fig. IS3I 
which are there indicated with the corresponding line style (see text 
for details). 
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FIG. S3: (Color online) Main panel: non-interacting magnetic sus- 
ceptibility xo (Lindhard function) for the quasi Id model. In the 
inset, Xo for the truly Id Hamiltonian (obtained when just the dom- 
inant hopping is taken in account, t'^±^ = 0.88). The arrows indi- 
cating peaks located at (fc, k) in the BZ are reproduced in Fig. IS4I 
These peaks indicate different spin fluctuations which may lead to 
electronic pairing. 

Figure|S2]shows the band structure for high symmetry lines 
of the BZ. The Fermi energy is located at Ep = 0.0. The 
two hole-pockets obtained are identical and nearly square (see 
Fig.|S4|i. their corners being slightly rounded due to the pres- 
ence of the finite t^^ = 0.05 hopping. For — 0.0 the 
hole-pockets are perfectly square and the two bands are de- 
generate along the F — M (S) line. 

Figure |S3] shows, in the main panel (solid red curve), the 
non-interacting magnetic susceptibility (Lindhard function) 
Xo obtained from the bands in Fig. [S2] The inset shows, as 
a reference, xo for the truly Id Hamiltonian (i.e., t^^ = 
0.0). A comparison of these two xo curves in Fig. |S3] with 
the one for the fully two-dimensional (2d) two-orbital model 
(Fig. [Ttb) in the main text) shows that the introduction of a 
small t^^' = 0.05 brings the Id model xa (green curve in the 
inset) qualitatively close to the 2d result. To see that, com- 
pare the solid (red) curve in the main pannel of Fig. IS 3 1 with 



the dashed (green) curve in Fig. [Ttb) of the main text. As 
will be described next, the extra peaks introduced in the S 
line (main panel of Fig. |S3t have a marked influence in the 
singlet pairing gap functions. Indeed, Fig. |S4ta) shows the 
dominant gap function (with symmetry Big) for n = 0.5 and 
J/U = 0.3. The four vectors connecting local maxima (with 
opposite signs) of the gap function are exactly the same that 
locate the four (fc, fc) peaks in xo in the main panel of Fig. |S3] 
This, once again, shows the strong connection between spin 
fluctuations and electron pairing. Panel (b) shows the sub- 
dominant gap function, with symmetry . Their eigenval- 
ues are the same up to the third decimal place. It is easy to 
see that the vectors displayed in the i?2g gap function [panel 
(a)] apply identicaly to the subdominant in panel (b). It 
is also interesting to observe that a possible extra set of pair 
scattering processes, leading to change of sign in the gap 
function, are the ones connecting adjacent sides of the same 
hole-pocket. Two of them are indicated by double-headed ar- 
rows. However, these processes do not occur, as there are 
no peaks in xo that can provide spin fluctuations with these 
two wave vectors (see Fig. IS3b . This results in the pairing 
strengths of both gap functions being basically the same. 

The presence of a relatively large number of different pair- 
ing spin fluctuations, as implied in Fig. IS4I suggests that the 
RPA spin susceptibility should have competing peaks when U 
is close to the critical value. Figure|S5]indicates that this is in- 
deed the case. There we show xrpa for the same parameters 
as Fig. |S4] for three different values of Hubbard interaction 
U = 0.77 (dotted blue curve), 0.79 (dashed green curve), and 
0.8 (solid red curve). As a comparison, at the same filling 
71 = 0.5 and J/U = 0.3, for the full 2d model studied in 
the main text, the leftmost peak (Fig. |2jb), main text), at a 
comparable ratio U/Uc as the ones in Fig.[S5] is a few orders 
of magnitude above the other peaks. Similar results are seen 
for the other fillings and J/U values, indicating that there is 
mainly a single dominant pairing process in the 2d model in 
the main text. In the Id model it seems as if the different wave 
vector spin-fluctuations cooperate to produce pairing. 
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FIG. S5: (Color online) RPA spin susceptibility for the quasi Id two- 
orbital model Hamiltonian, n = 0.5 and J/U = 0.3. Three curves 
are shown for values of Hubbard U close to the critical value Uc ~ 
0.8686; U = 0.77, 0.79, and 0.8. A competition between two peaks 
can be clearly observed. They are indicated by the same type of 
arrows as the ones for the corresponding peaks in Fig |S3l 



